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Abstract 

We calculate and analyze Z and production in association with quark-anti- 
quark pair in fc^-factorization. Numerical calculations are performed using the Monte 
Carlo generator Cascade for proton proton collisions at LHC energy. We compare 
total and differential cross sections calculated in fey-factorization approach with total 
differential cross sections obtained in LO and NLO calculations in collinear factoriza- 
tion approach. We provide strong evidence that some of the effects of the NLO and 
even higher order collinear calculation are already included in the LO /cr-factorization 
calculation. 



1 Introduction 

In the following years new discoveries are expected at the LHC concerning physics within 
the Standard Model and beyond it. The discovery of the Higgs boson and exclusion or 
affirmation of possible extensions or alternatives to the Standard Model will be of special 
interest. To be able to measure the proposed signals of processes which open the access to 
new physics a very good understanding of the detectors and their responses to produced 
particles will be needed. An accurate calibration of particle detectors could be achieved 
by using processes with well known cross sections in which particles with well known 
properties are produced. A calibration of LHC detectors using W or Z signals is proposed 
in several publications [1]. Moreover, the W or Z production is important because it 
plays a significant role in background processes connected to Higgs production. Another 
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experimental motivation is provided by the possibility to measure the luminosity via Z 
boson production [2]. 

At the Tevatron collider W/Z production takes place at a typical x = \[M^Js ~ 0.04 
and hence is dominated by scattering of quarks. Because of the much higher energy, proton 
scattering at LHC will allow smaller proton energy fractions and will be dominated by gluon 
scattering. 

The W mass provides a hard scale and allows a perturbative calculation of the hard ma- 
trix element. The resummation of large logarithms of the form [a s ln(/i 2 /AQ CD )] n (where 
/i 2 ~ M^, jj 2 ^> Aq CD ) can be performed in the framework of the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) equation [3], leading to the collinear factorization into 
conventional parton densities and a hard scattering matrix element. While in the conven- 
tional collinear approach the longitudinal momentum fraction is considered to be dominant, 
such that the transverse momenta of the partons can be neglected as well as their virtual- 
ities, at small x the transverse momenta entering the hard matrix element should become 
relevant. 

At the LHC the larger center of mass energy allows W/Z production at even smaller x 
such that the production of particles will be dominated by gluon-gluon fusion. Moreover, 
in this situation we have to deal with two different large scales (s ^> jj 2 ^> Aq CD ) and 
logarithms of the form [a s \n(l/x)] n arise which have to be resummed. This is realized 
by the leading logarithmic (LL) Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation [4] or 
the Ciafaloni-Catani-Fiorani-Marchesini (CCFM) evolution equation [5] which additionally 
resums terms of the form [a s ln(/i 2 /AQ CD )] n and [a s ln(/i 2 /AQ CD ) ln(l/a;)] n . Just as for 
DGLAP, it is possible to factorize the cross section into a convolution of process-dependent 
hard matrix elements with universal parton distributions. But as the virtualities and 
transverse momenta are no longer ordered (as it is the case in DGLAP evolution), the 
matrix elements have to be taken off-shell, and the convolution has to be made also over 
transverse momenta with the so-called unintegrated parton densities. This factorization 
scheme is called kx -factorization [6, 7] or semi-hard approach [8] and will be used in this 
work. 

There is also the notion of transverse momentum dependent (TMD) parton distributions 
[9]. But although in these approaches the transverse momentum of the parton is taken into 
account as well, this is only the case on the side of the parton density. The matrix element 
is calculated with incoming on-shell partons, and transversal momenta of the incoming 
partons are neglected. It has been shown [10] that factorization within this approach is 
violated beyond NLO. In case of the /^-factorization approach used in this work this is 
also expected. Indeed, it is well known that in the BFKL approach beyond NLO multiple 
gluon exchange in the t-channel has to be taken into account. 

In this paper we calculate and analyze Z and W production associated with two quark 
jets provided by gluon-gluon fusion in /cy-factorization. We assume quasi-multi-Regge- 
kinematics (QMRK) where the cluster of W/Z and the two quarks is well separated in 
rapidity from the proton remnants while the kinematics within that cluster is considered 
without any further assumption. In particular, we take into account the mass of the 
quarks. In this kinematic regime a gauge independent off-shell matrix element can be 
extracted due to high energy factorization. A similar calculation has been done in [11], 
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Figure 2.1: Labeling and flow of momenta of the process pp q (W/Z) qX. 



where the authors calculated photon (instead of Z/W) production in the same framework. 
We calculated the matrix element independently and extended it to massive gauge bosons. 
In our work on massive gauge bosons production we especially focus on the predictions for 
LHC and compare with a collinear factorization based calculation. 

When this paper was in preparation, we learned about another group [12] working on 
this process as well using the same theoretical approach, but laying more emphasis on 
confronting the theoretical predictions with experimental data and examining the role of 
quark contributions. 

The paper is organized in the following way: In section 2 we describe notation, kine- 
matics of the process and the calculation of the matrix element. In section 3 we present nu- 
merical results obtained from a calculation using the Monte Carlo generator Cascade [13], 
where the matrix element squared was implemented. In section 4 we summarize the results 
and offer conclusions. 



2 Kinematics of Z/W production and calculation of the 
hard matrix element 

We label the 4-momenta of incoming hadrons with masses rriA and mg by p' A and p' B , 
respectively. In the center of mass system they can be expressed in terms of invariant light 
like vectors pa and ps 

2 2 

Pa =Pa + Pb, Pb =Pb + Pa- (2.1) 

s s 

In the case of protons at the LHC we have m 2 A = m 2 B = m 2 which satisfies the relation 

2 

Y < 1. Therefore, we can neglect the masses in Eqs. (2.1) and use Pa,b instead of p' AB . 
It is convenient to use Sudakov decomposition for all momenta present in the calculation 
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(see also Fig. 2.1 1 ) by decomposing them into components proportional to pa and pb, and 
a remainder perpendicular to both of them 

ki = aip A + PiPB + k±, (2.2) 

where i G {1, 2, W(Z)} for outgoing particles, and 

qi =ap A + (3 qi PB + qi±, q2 =a q2 PA + Ppb + <?2_l (2.3) 

for the gluons entering the hard matrix element. It is also convenient to introduce Euclidean 
two dimensional vectors ki and cfj which satisfy the relations k 2 = —k 2 ± > and q 2 = 

-q] ± >0. 

In QMRK we have 

<fi = -tf = ti, (2.4) 
/3»a w , q l = -ql=t 2) (2.5) 

aA =™i±M. i (2.6) 
s 

where % G {1,2, W(Z)}, and m ; are the corresponding masses of outgoing particles. The 
invariants t\ and t 2 describe the momentum transfer between the cluster formed by the 
quarks and the W (Z) boson on one hand and the incoming protons on the other hand. 
Due to the strong ordering in a and (5 one can neglect terms proportional to f3 qi and a q2 
in the calculation. 

It is useful to introduce a set of Mandelstam variables describing the system 

s =(gi + g 2 ) 2 = a(3s - (ft + q 2 ) 2 , (2.7a) 

§! =(h + k w ) 2 , s 2 ={k 2 + k w ) 2 , (2.7b) 

ii =(Qi ~ k i)\ i 2 =(q 2 - k 2 ) 2 , (2.7c) 

ui =(9i - h) 2 , u 2 ={q 2 - h) 2 , (2.7d) 

related by 

ui + i 2 + s = ti+t 2 + ml + s 1} u 2 + i 1 + s = t 1 +t2 + ml + s 2 . (2.8) 
It is convenient to introduce transverse masses defined by 



m i± =\Jm 2 + k 2 , m q± =a/s + (gl + q 2 ) 2 , (2.9) 

and longitudinal momentum fractions of the produced particles Xi — — . Combining these 
relations with Eqs. (2.6, 2.7) one finds that - in the end - the matrix element of W or Z 
production associated with a quark-antiquark pair can be expressed in terms of independent 
Mandelstam variables defined in Eqs. (2.7), transverse masses and variables x lj2> w(z)- 



1 These and the following diagrams were drawn in JaxoDraw [14]. 
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In the fcr-factorization formalism the hadronic and partonic cross section are related 
as follows: 

da(pp^q(W/Z)qX) = J ^ J dqt J ^A{a,qt,f) 

X jjj d ^l ^A(P,<f 2 ,fz 2 )d&(g*g* - q(W/Z)q), (2.10) 

where A is the unintegrated gluon density in a proton and 01,2 is the angle of q\ t 2 with 
respect to some fixed axis in the azimuthal plane. The argument /i 2 of unintegrated gluon 
densities is the factorization scale. The partonic cross section is denoted by da. 

Since the incoming gluons of the matrix element entering this partonic cross section 
are off-shell, the calculation differs from that of a hard matrix element in the collinear 
approach significantly. To guarantee gauge invariance, the process with off-shell incoming 
particles has to be embedded into the scattering of on-shell particles. The extracted off- 
shell matrix element is of course independent of the specific choice of the particles in which 
the scattering process is embedded. Therefore, we replace the protons by quarks for the 
calculation of the hard matrix element. All diagrams for the discussed process are shown 
in Fig. 2.2. 




Figure 2.2: Full set of diagrams contributing to W/Z production via off-shell gluon-gluon 
fusion. 

The first two rows of Fig. 2.2 include also non-factorizing ('non-resonant') diagrams 
which factorize only in the sum. To make this factorization apparent already at this level, 
one can sum up the different diagrams of one gluon production in quark-quark scattering 
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leading to one effective diagram with an effective vertex (see Fig. 2.3). By working in 
Feynman gauge one obtains the well known Lipatov vertex [15]: 



2p Ar p Ba ( 2h + m 2 q± 2t 2 + m 2 q± \ 

' «&) = {—fa— p A ~ ^^ PB ~ iqi± ~ ^ ) ■ (2 - U) 

It can be shown that this vertex obeys the Ward identity. By this procedure, the first two 
rows of Fig. 2.2 are each replaced by just one diagram. 




Figure 2.3: Diagrams contributing to the Lipatov vertex. 



Strong ordering of Mandelstam variables s and tip allows us to make a simplification 
of the coupling of gluons to incoming quarks. By neglecting the exchanged momentum 
in the vertex, we get an eikonal vertex which does not depend on the spin of the particle 
coupled to gluon and preserves its spin. In detail, it reads 

-iu^pA-qi^uiKpA) — > -2ip» A 5 KM . (2.12) 

With the help of Eq. (2.12) it is possible to remove the external quark lines and attach 
so-called 'non-sense' polarizations to the incoming gluons: 

V 



e" e" (2 13) 



Instead of Feynman gauge, one can choose an appropriate axial gauge [6] n ■ A = with 
the gauge vector 

n v = ap\ + bp B with a, b e C. (2.14) 

The contraction of the eikonal coupling (2.12) with the gluon polarization tensor in this 
gauge 

<>te) = -a. + 'M^J^ _ „ 2 ^_ (2 . 15) 

then reads 

P>^) ~> P^fo) =^ (2.16) 
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In such a physical gauge the 'non-resonant' diagrams vanish since the direct connection 
of two eikonal couplings gives p^d^Jps = (in other words: the Lipatov vertex is to be 
replaced by the usual three gluon vertex). 

In the case of heavy quark production the polarization sum for the s-channel gluon 
reduces to its Feynman gauge analogue —g^v due to the heavy flavor current conservation. 
The same simplification takes place in our calculation. Nevertheless, we have to stress 
that in general the polarization sum stays in its complex form. Of course, both ways to 
calculate the matrix element are equivalent due to gauge invariance. 

The sum over the physical polarizations 77 of the W boson reads 

V 6^(77, k w y»( v , k w ) = -gr + (2.17) 

^ m w 
It is equivalent to replace the polarization sum by 

J2^(v,k w )e* u (v,k w ) - -<T, (2-18) 
v 

and to add also the contribution of the Goldstone boson emission diagrams, where the 
W boson is replaced by a Goldstone boson with mass mw- This is in analogy of using 
the Feynman-t'Hooft gauge instead of the unitary gauge. We have calculated the squared 
matrix element in both ways as a crosscheck. 

Expressions for the single diagrams in Fig. 2.2 - where the first diagrams are already 
combined using the Lipatov vertex of Eq. (2.11) - are listed here (the hat marks contraction 
with Dirac-matrices): 

KA ab ■ 2j7 -\.b .allV^l,^) ~fa - k W + mi 5 

M lliv = -tg w g s K w/z u[t ,t ] 2 e K ~ a <?7 K 

S S 2 — Tfl-y 

KA ab ■ 2ts + kw + m 2[b , a] t^(q 1 ,q 2 ) 
M 2fMl/ = -ig w g s K w/z ue(v q - a q Y) : 2 — [t ,t ]-^-; v, 

KA ab ■ 2t, - + a ki-qi + m lb -fa-kw + rrii 5 
Mf^ = -ig w g s K w/z ut 7 M ^ 2 — 1 lv : 2 e(v q - a q Y)v, 

t\ — fflj S 2 771 j 



k a (ib ■ 2 7v ^ki + k w + m 2 q 2 -fa + m 2±h 

M 4 = -tg w g s K w/z ue{v q - a q ^ ) — — t 7 M — — t 7^, 

si — m 2 t 2 — 777,2 

KA ab - It, -,a fa - q x + m l 5 q 2 - fa + m 2 b 
M blxv = -ig w g s K w/z ut — e(v q - )— — t 7^, 

KA ab ■ 2t, -~i 5\fa + k w + m 2 b qi-k 2 + m 2 
M 6fU/ = -ig w g s K w/z ut{v q - a ql ) - — 2 — t lv — -5— t 7^, 

KA ab ■ 2r, -,b fa-fc + mi.a ~fa - k W + 7?7i 5 

M 7flu = -ig w g s K w/z ut lv — -3— t 7m , — 2 e(v q - a ql )v, 

Ui 2 lll-y h 2 III^ 

K,ab ■ 2r/ -a fa-q 2 + mi ^q\-fa + m 2 
M a l„ = -ig w g 2 s K w/z ut b lv e(v q - a ql 5 )\ -5-^, 
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(2.19) 



with the short hand notations u = u(X, k\), v = v(X f , k 2 ), e = e(r), kw), and where rj, A 
and A' label the helicity/ spins of the corresponding particles. Color factors are represented 
by Gell-Mann matrices t a , t b . The factors v q , a q and K w / Z encode the W and Z coupling. 
For W boson we have v q — a q — 1 and Kw = VudT^, where V u d is the corresponding 
element of Cabibbo-Kobayashi-Maskawa matrix. For Z we have a u — |, v u — | — | sin 2 9 W 
and da = — ^, Vd — — | + | sin 2 W and = 2c ole w ' w here 6*^ is the Weinberg angle. In 
the latter case mi equals m 2 , and miy is replaced by mz- 

If we make use of the Eq. (2.18) to replace the polarization sum, one has to add diagrams 
and corresponding amplitudes with Goldstone bosons with couplings 

. {1712-171! m 1 +m 2 5 \ , . 

-ig w K w/z [ v q ag'j . (2.20) 

V mw/z rnw/z ' 

Finally, the square of the amplitude averaged over initial helicities and colors of gluons 
and summed over spins/ helicities and colors of final particles can be written as 



4 (N* - l) 2 1 1 4 (TV? - 1) 

V c ' \ c / x,\',T),a,b 



color ' 



c gi c <?2- / 1 «a 



i ab 

1=1 



(2.21) 



By evaluating the traces over the products of Gell-Mann color matrices, one encounters 
two possible cases of color factors 

i yv 2 — 1 1 ( N 2 — 1 ) 2 

Tr{t a t b t a t b } = - - c — , Tr{t a tHH a } =- 1 c — } , (2.22) 



c 



where N c = 3 is the number of colors. 

Finally, the expression for the partonic off-shell cross section appearing in Eq. (2.10) 
to calculate the hadronic cross section is 

da(g*g* - q (W/Z) q) =(2tt) 4 (5^ (q, + q 2 - h - k 2 - k w/z ) x 

d% (2.23) 



2a(3s t ± t 2 4(N 2 -iy l 1 H (27r) 3 2E(A; l ) ' 
r 1 2 v c 7 ie{i,2,w(z)} v y v 1 



The origin of the specific form of the flux factor and prefactor a f lt * is formulated in 
[6,7]. We summarize the most relevant aspects here. An important feature of the whole 
calculation is that it is possible to recover the result obtained in collinear factorization by 
neglecting the transverse momenta of the gluons when they enter the hard matrix element 
and instead integrate over them only in the gluon densities. Due to factorization it is pos- 
sible to keep this connection not only for the full cross section, but also for gluon densities 
and hard matrix element separately as well, provided that the explicit manifestations of 
the factorization formulae are phrased. 

The key point is the observation that 

2^f^\ =-^i = (2^f^\ . (2.24) 
qj ± /*i \ q\ L /<t>2 V ' 
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As shown in Eqs. (2.14-2.16), in an appropriate gauge the polarization sum 2va ^ Bv can be 
replaced by 2qiJ a ^ ± " ■ Since in this gauge one has to deal with exactly the same diagrams 
as in the on-shell calculation, by dressing the off-shell matrix element squared with the 
prefactor s and performing the averaging over azimuthal angles of the 'incoming' 
gluons, followed by taking the limit ti,t 2 — > 0, one gets the collinear limit of the matrix 
element squared. The flux factor for off-shell gluons is defined as for on-shell gluons with 
2^jgj- As the matrix element is gauge invariant, this connection remains valid when one 
performs the current calculation in a different gauge. 

Due to the off-shellness of the incoming gluons and the three particle final state the 
final result of the matrix element squared is rather lengthy. For that reason, we calculated 
it independently and in different ways. One calculation followed directly the derivation 
above using Feynman gauge for the gluons, and has been performed using M ATHEMATICA. 
A second calculation written in Form [16,17] used an axial gauge as described above such 
that the Lipatov vertices in (2.19) are to be replaced by standard three-gluon- vertices. 
Moreover this second method used the method of orthogonal amplitudes, described in [18], 
which affects the fermionic part of the matrix element and with which one is able to treat 
the matrix element squared in a more compact way. 2 

For this second method a few technical details are elaborated in the remainder of this 
section. The method of orthogonal amplitudes is based on expressing a generic amplitude 
M. (with one quark line) in terms of a set of four independent operators Oi, % G {1, ..,4}, 

which satisfy orthogonality relations Tr{Oj(A;2 ~ m 2 )Oj(^i + mi)} = \\Oi\\ 2 5ij for any 
possible i and j, where ||0,|| is the "norm" of the operator Oj. The projection of Ai by an 
operator Oi is performed in the following way 

M = -^^Mv{X',k 2 )Oia{\M)- (2.25) 

11^*11 A, A' 

The matrix element squared then has the following form 

\M\ 2 = \Mf ■ (2.26) 

A, A' i 

In our case the matrix element consists of up to five Dirac-matrices (neglecting 7 5 ), after 
squaring one has to evaluate traces of up to twelve of them. In contrast the method of 
orthogonal amplitudes leads only to traces of up to eight Dirac-matrices. 

If one wants to consider also the Z or W ± coupling in the Feynman diagram, one 
encounters a technical problem connected with the appearance of the Dirac-matrix 7 5 in 
the expression for the amplitude, leading to terms which include Levi-Civita tensors which 
later cancel. To avoid this complication, one can split the expression for the amplitude 
into two parts, one which does not include 7 5 and the other one which does (to separate 
the vector and axial part of the Z or W boson coupling). For the part with 7 5 one uses a 
base of operators 0«7 5 . It is easy to check that they satisfy the same orthogonality relation 



2 We also have cross-checked numerically our results for the case of a produced photon instead of a W/Z 
boson with those of the authors of [11] whose cooperation we gratefully acknowledge. 
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like the operators 0{. One also easily see that projections of amplitudes in which 7 5 occurs 
do not contain terms with Levi-Civita tensors. In doing so, we extend the method of 
orthogonal amplitudes in a natural way. 

Another complication comes from the presence of color factors in the expressions which 
are not numbers but matrices. To treat the projections as numbers, it is necessary to 
separate the Feynman diagrams into three groups according to different color factors, 
namely 

C ab = t a t b _ t b f a^ 

Cf = t a t\ (2.27) 
Cf = t b t a , 

which form a vector C ab = (Cf , Cf, Cf) (components of C ab are color factors of M^ t2)fJH/ , 
•^f3-5W an d -M?6-8W correspondingly). One can then build a corresponding vector 
containing the sums of Feynman diagrams without the color factors T = (J-i, JF 2 , JF 3 ) such 
that 

M ab = (C-ob^T^ ( 2 _ 28 ) 

The Lorentz indices have been dropped for simplicity. Using the matrix 

dj = Tr{CfCf}, (2.29) 

the expression for the square of the matrix element takes the form 

\M\ 2 = 7*CF, (2.30) 

where combinations of Ti and T* are calculated using the projection method introduced 
in Eqs. (2.25, 2.26). For the final simplification we have diagonalized the matrix C. After 
diagonalization of the matrix C only two diagonal elements remain nonzero. This is ex- 
pected because the quarks in the final state, in this process, can occur only in two possible 
color states. 



3 Numerical studies 

The last missing pieces needed to calculate the hadronic cross section using Eq. (2.10), are 
the unintegrated gluon densities. As mentioned in the introduction, there are two equations 
suited to describe the evolution of an unintegrated gluon density, namely BFKL [4] and 
CCFM [5], respectively. Both have been shown to agree on the leading logarithms in small 
x [19], but the CCFM evolution is valid in the domain of larger x as well and, moreover, 
matches in this region with DGLAP. Therefore, we base our numerical studies on an 
unintegrated gluon density obeying the CCFM equation, which has been implemented in 
the Monte Carlo generator Cascade [13]. We also investigate how the results change 
when using uPDFs generated by a different procedure known as KMR [20]. 

For this purpose, we implemented the matrix element squared as described above into 
Cascade. This implementation will be available in the next version of Cascade. 
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We have used the unintegrated parton distribution function (uPDF) CCFM 2003 set 3 
for the numerical calculation. 

To investigate the calculated matrix element as accurately as possible, we neglect in 
this first study the effect of hadronization of the final state. We study in detail rapidity 
and transverse momentum distributions of the produced gauge boson, quark and antiquark 
which (if one assumes that quarks approximately determine jets) are the most important 
observables in the experiment. 

Furthermore, we compare the /^-factorization approach to the collinear one. For this 
purpose, we compare the distributions obtained by our transverse momenta dependent ma- 
trix element with distributions obtained from the Monte Carlo generator Mcfm [21] which 
provides a calculation of the same process in the collinear limit. In that case the transverse 
momenta coming from the evolution are neglected. We also investigate in Sec. 3.3 how 
the variation of unintegrated parton densities affect the azimuthal angle and transverse 
momenta distributions. 

As an artefact of the perturbative calculation, the results depend on the renormalization 
scale and the factorization scale /if. In the CCFM formalism the hardest scale is set by 
the emission angle of the hardest subcollision. It can be expressed in terms of the energy 
of the subcollision /if — a/s + (gl + q^) 2 - For the comparison with collinear factorization 
calculations we have used as renormalization scale /jLr = mz in /^-factorization calculation 
and in collinear calculation as well. We have also investigated other possible choices (see 
subsection 3.3). 

3.1 Comparison with LO collinear calculation 

Our calculation of the hard matrix elements includes W ± and Z production in association 
with all possible quark-antiquark channels in gluon gluon fusion. Since the basic structure 
of all these matrix elements is very similar, we present results only for the typical case 
of Zbb production at LHC energies of y/s = 14TeV. The mass of the 6-quark used is 
rat, = 4.62 GeV. For the collinear factorization calculations we use the parton densities 
CTEQ6L1 [22]. 

The total cross sections are comparable in magnitude, though they differ considerably: 
0.406 nb in /c^-factorization and 0.748 nb in collinear factorization. The difference of total 
cross sections stems from the different behavior at low transversal momenta of final state 
particles (discussed later in this section) where contributions from transversal momenta of 
the initial state gluons play a significant role. It can be seen that by applying a cut on 
the transversal momentum of the Z boson pz± > 50 GeV the difference of the total cross 
sections becomes smaller. With this additional cut one obtains cross sections of 0.118 nb 
in /^-factorization and 0.141 nb in LO collinear calculation. 

The total cross sections for other final states of interest are given in Tab. 3.1. 

The transverse momentum and rapidity distributions of the vector boson are shown 
in Fig. 3.1 and 3.2, respectively. The comparison of the /^-factorization approach to the 
collinear shows that they agree in transversal momentum distributions of Z at high values 
of this quantity. This is no surprise, since at high pz± the contribution from initial state 
gluon transverse momenta is expected to become small. 
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P Zi [GeV] 

Figure 3.1: Transverse momentum distributions of the produced Z gauge bosons. Calcu- 
lation with massive 6-quarks. Both calculations are in LO of perturbation series. 




Figure 3.2: Rapidity distribution of the produced Z gauge bosons. Calculation with mas- 
sive 6-quarks. Both calculations are in LO of perturbation series. 
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final state 


Zee 


Zbb 


Ztt 


W + sc, W'cs 


otot [nb] 


0.430 


0.406 


0.525- 10- 3 


1.92 



Table 3.1: Total cross sections for different final states, calculated in fc^-factorization using 
Cascade. 




The rapidity distributions of the Z show a similar behavior, except for the overall 
normalization (Fig. 3.2). 

To elaborate the difference between hp- and collinear factorization, we investigate more 
exclusive observables, like the cross section differential in rapidity distance between quark 
and antiquark (Fig. 3.3). Both calculations show a two peak structure with a minimum at 
zero rapidity, but the /^-factorization result has a considerably shallower minimum. The 
minimum in the case of the collinear calculation gets shallower - bringing together both 
calculations - when one again applies a cut on pz± > 50 GeV as one can see in Fig. 3.4. 

In the distribution of the azimuthal angular distance of Z and max^i,^^) (Fig. 3.5) 
we observe that the region from to 7r/2 is forbidden within the collinear calculation due to 
momentum conservation, which is not the case for fey-factorization. This is caused by the 
contribution from initial state gluon transversal momentum which allows the transversal 
momenta of Z, b and b to be unbalanced. A larger spread of possible configurations causes 
that the distribution in the /cr-factorization calculation flattens. 
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Figure 3.4: Distributions of the rapidity distance between quark and antiquark. Calculation 
with massive 6-quarks. A cut on pz± > 50 GeV has been applied. 




Zhb 



Figure 3.5: Distributions of the distance in azimuthal angle of Z and highest p± quark or 
antiquark. Calculation with massive 6-quarks. Both calculations are in LO of perturbation 



series. 
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3.2 Comparison with NLO collinear calculation 

In collinear factorization the physical effect of the intrinsic transverse momenta of the initial 
gluons can not be described until higher order corrections are taken into account. Then 
additional real emissions lead to off-shell gluons and their transverse momenta. Therefore, 
the significant differences between a calculation in the collinear factorization framework and 
/cT-factorization framework shown in the previous section encourage us to compare our LO 
calculation in /cr-factorization with a NLO collinear calculation, since CCFM evolution 
includes the high-energy part of the NLO collinear corrections. Since there are two off- 
shell initial gluons in a /cr-factorized pp-collision, one could even call for a higher order 
collinear calculation to compare 3 . 

To compare with a collinear NLO calculation, we use again the Monte Carlo generator 
Mcfm. This Monte Carlo generator provides the process gg — > Zbb at NLO only in the 
massless quark limit. To avoid divergences, additional cuts are applied on transversal 
momenta of quarks, on the invariant mass of the bb pair, and on transversal momenta 
of a gluon which is produced in diagrams of real NLO corrections. Transversal momenta 
of produced quark, antiquark and gluon have to satisfy the condition p± > 4.62GeV 
(corresponding to the mass of the 6-quark). These cuts on quark (antiquark) momenta are 
automatically applied in Mcfm when one is performing a calculation involving massless 
quarks (antiquarks). We choose the parton density functions set CTEQ6M [22]. The same 
cuts on transversal momenta of quark and antiquark are then applied in Cascade as well. 

For the total cross sections, we obtain in the NLO collinear factorization calculation 
1.04 nb, and in the /c^-factorization calculation 0.429 nb. The difference of the total cross 
sections in /^-factorization calculation and the NLO calculation in collinear factorization 
is of the same origin as the difference between the total cross sections in section 3.1 where 
comparison of /cr-factorization calculation and NLO calculation in collinear factorization 
is discussed. This is again illustrated by a cut on pz± > 50 GeV diminishing the difference 
between the cross sections (0.125 nb for the /^-factorization calculation and 0.165 nb for 
the NLO calculation in collinear factorization). 

The result for the cross sections differential in the transversal momentum of Z can be 
seen in Fig. 3.6. The cross section changes especially at small pz± (see Fig. 3.7) from LO 
to NLO calculation, and the difference between collinear calculation and /cr-factorization 
calculation becomes more pronounced. We observe that the maximum of the distribution in 
the NLO calculation (Mcfm) stays approximately at same value of transversal momenta 
and the shape of the peak is very different from the one we obtain in /^-factorization. 
Nevertheless, the pz± distributions match at very high pz± (O(10 2 GeV)). 

The rapidity distribution of the Z (Fig. 3.8) shows no major difference in shape in 
//^-factorization approach, LO and NLO collinear factorization approach. 

3 Although we argue that already the LO /cr-factorization calculation includes in some sense higher order 
corrections, one might ask for an extension to NLO. So far fcr-factorization based on CCFM evolution 
has been formulated only at LO. On the other hand, since the BFKL equation has been calculated at 
NLO accuracy [23], in the small x regime fcr-factorization can be formulated at NLO accuracy as well 
[24]. Nevertheless, an implementation into a Monte Carlo generator is still outstanding. Moreover, the 
calculation of an off-shell 2 — ► 3 process at one loop order is far beyond the scope of this work. 
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Figure 3.6: Comparison of cross sections differential in transverse momentum of the pro- 
duced Z gauge boson. Calculation with massless 6-quarks. The applied cuts are described 
in the text. 
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Figure 3.7: Comparison of cross sections differential in transverse momentum of the pro- 
duced Z gauge boson (linear scale). Calculation with massless 6-quarks. The applied cuts 
are described in the text. 
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Figure 3.8: Comparison of cross sections differential in rapidity of the produced Z gauge 
boson (logarithmic scale). Calculation with massless 6-quarks. The applied cuts are de- 
scribed in the text. 



We consider the cross section differential in the total transversal momentum of the Zbb 
system Pzblx m Fig. 3.9. In the NLO collinear calculation a non-zero p Zb b± i s generated by 
the emission of an additional gluon, while at LO it is always balanced to zero. At low p Zb b± 
we see the consequence of the cut on the transverse momenta of the outgoing particles in 
Mcfm (a small gap between GeV and 4.62 GeV in p Zb i± histogram). Since there are no 
parton showers or soft gluon re-summation [25] included in the Mcfm NLO calculation, 
one observes a steep rise of the cross section towards zero transverse momentum because 
the matrix element diverges when approaching p Zbb - ± -^>0 GeV . On the other hand, uPDFs 
include corrections similar to parton shower effects, treated consistently, which causes the 
turnover in the cross section of the fc^-factorization calculation. Here, the entire transversal 
momentum of the Zbb system stems from the transversal momenta of initial state gluons. 
We expect that resummation effects at low values of p Z bi± would tame the growth of the 
cross section in collinear factorization and would decrease the difference to fc-r-factorization. 
Interestingly, there is a difference not only at low values of p Zbb - ± , but also at high values 
of 

Pzblx- The differential cross sections at high p Zb i± have a similar slope, but differ by a 
factor of ~ 3. This is contrary to the behavior of distributions of p z ± in Figs. 3.6 and 3.1 
where at large values of p z ± the differential cross sections overlap. For this difference at 
large p Zb b± further calculations have to reveal the exact effect of higher order corrections 
in collinear factorization, keeping in mind that the NLO for this obersvable de facto is the 
first non trivial order. 

The cross section differential in the difference of azimuthal angles of Z and b or b quark 
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Figure 3.9: Comparison of cross sections differential in the p± of the system Zbb. Calcula- 
tion with massless 6-quarks. The applied cuts are described in the text. 



with higher transversal momentum - A<pzhb - is shown in Fig. 3.10. Going from LO to 
NLO, the collinear calculation reveals a broader distribution like in the /c^-factorization 
case. Nevertheless, the fcr-factorization result shows a more homogeneous spread of the 
azimuthal angle distance. This difference origins partly in the difference of the transversal 
momentum distributions at low values (see Fig. 3.7). A cut on low values (pz± > 50 GeV) 
of the transversal momentum of the Z boson results in steeper A(p Z hb distributions as 
shown in Fig. 3.11. Still, the fcy-factorization result is flatter than the NLO collinear 
factorization calculation giving an indication that there is a contribution from the total 
transversal momentum of the Zbb system generated by both uPDFs. 

3.3 Variation of the CASCADE results on uPDF and renormaliza- 
tion scale 

To estimate the uncertainty coming from the different choices of uPDF sets, we calculate 
the cross sections differential in either the transverse momentum of the Z boson or Aifizhb 
(distance in polar angle between Z and max(p;,x,p5jj) using different sets of uPDFs, namely 
CCFM J2003 set 1, 2, 3 [26] and CCFM set AO [27], which are all obtained from fits to 
HERA F 2 data [28]. In addition we use the unintegrated parton density by [20], referred 
to as KMR. The resulting plots are shown in Figs. 3.12 and 3.13. We do not show the 
distributions for set 1, because they are very close to distribution for the set 3, to keep the 
plot clear. 

The total cross sections obtained for different uPDFs can be seen in Tab. 3.2. The 
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Figure 3.10: Comparison of cross sections differential in distance in azimuthal angle of Z 
and higher p± b/b. Calculation with massless 6-quarks. The applied cuts are described in 
the text. 




Figure 3.11: Comparison of cross sections differential in distance in azimuthal angle of Z 
and higher p± b/b. Calculation with massless 6-quarks. An additional cut on pz± > 50 GeV 
has been applied. 
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Total cross section [nb] 


PPFM 7900^ «p+ 1 

vy^r 1V1 JZUUO ocl 1 


u.ouy 


CCFM J2003 set 2 


0.147 


CCFM J2003 set 3 


0.406 


CCFM set BO 


0.277 


CCFM set AO 


0.378 


KMR 


0.190 



Table 3.2: Total cross sections of the process pp — > Zbb+X for different sets of unintegrated 
parton distribution functions. 





Total cross section [nb] 


m z 


0.406 


2m z 


0.392 


\m z 


0.607 


V m l+P 2 z± 


0.467 


2y/m 2 z +p 2 z± 


0.381 


W™z+Pz± 


0.585 



Table 3.3: Total cross sections for different renormalization scale //. 

total cross section varies for these different uPDFs about 45%, while the shape of the 
distributions is hardly effected except of the KMR. KMR set uses completely different 
evolution equations and a deviation is not surprising. 

As a last point to discuss, we turn to the scale dependence. As already mentioned in 
the beginning of section 3 the factorization scale is fixed by the emission angle of the hard 
subprocess. However, there is still freedom in choice of the renormalization scale which 
should be of order of the typical scale of the hard subprocess. 

We consider two possible choices: the constant renormalization scale /xi = m z and the 
scale /i 2 = \/ m l + P%±i which are varied by factor of 2, so \x has values 2// 1; ^/ix and 2/j, 2 , 
\H2- The results for the p z ± and the A<j>zhb distribution can be seen in Figs. 3.14 and 
3.15, respectively. The values of the cross section for individual choices of the scale are 
summarized in Tab. 3.3. One can see that a running as does not affect the shape of the 
distributions, but only the total cross section. 

4 Summary and Conclusions 

In this paper we have calculated the matrix element for the process g*g* — > W/Zqiqj, 
taking into account the virtuality and transversal momenta of the initial gluons in the kr- 
factorization formalism. We have implemented the matrix element squared in the Monte 
Carlo generator Cascade and have calculated the total and differential cross sections of 
this process in proton proton collisions for the LHC at energy of = 14TeV. We have 
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Figure 3.12: Transverse momentum distributions of produced Z gauge boson calculated in 
Cascade using massive quarks. Cases with different uPDFs compared. 




Figure 3.13: Comparison of cross sections differential in distance in azimuthal angle of Z 
and higher p± b/b, using massive quarks. Cases with different uPDFs compared. 
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Figure 3.14: Transverse momentum distributions of produced Z gauge boson calculated in 
Cascade using massive quarks. Cases with different renormalization scales hr compared. 




Figure 3.15: Transverse momentum distributions of produced Z gauge boson calculated in 
Cascade using massive quarks. Cases with different renormalization scales hr compared. 
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compared our results with results obtained in collinear factorization (using Mcfm). The 
total cross sections differ by a factor of ~ 2. There are differences in distributions which 
are sensitive to compensation of transversal momenta of particles in the final state coming 
from rather fundamental differences between the two approaches. 

We found the most significant differences in the cross section differential in the az- 
imuthal angle between the Z boson and higher pt quark or antiquark - Acftzhb- While 
for a LO calculation in collinear factorization a region of values of Acpzhb is kinematically 
forbidden, in /^-factorization the whole range of A(ft Z hb is allowed. This is because of 
neglecting the contribution of transversal momenta of initial state gluons in calculation of 
matrix element in collinear factorization. The NLO collinear calculation (where transversal 
momentum is generated by real corrections) shows already the same qualitative behavior 
as the /cT-factorization calculation. However, there remains a difference in the shape of the 
distribution of A(pzhb compared to the /cT-factorization calculation We also compared cross 
sections differential in the transversal momentum of the Zbb system - Pzbb±- I n collinear 
factorization and lowest order perturbation theory (a|), the observable PzbB± is exactly 
zero. For a non-zero contribution in collinear factorization higher order corrections are 
needed. The //^-factorization gives non-zero contribution already at a| order. We have 
compared cross sections differential in p Z bb± calculated in NLO in collinear calculation and 
LO in /c^-factorization. The distributions have different shape at low values of PzbB±- 

At 

high Pzbl± the slopes are similar but differ in absolut size. 

We have calculated the cross sections differential in the transversal momentum of the 
produced boson. The maximum of the distribution in the /^-factorization calculation is 
at higher transversal momenta compared to the collinear one. This shows the sensitivity 
of this distribution on parton evolution model and treatment of kinematics. 

We conclude that some of the effects of NLO and even higher order collinear calculation 
are already included in the LO /cT-factorization calculation. 
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